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ELECTRON AVALANCHE STRUCTURE DETERMINED BY RANDOM WALK THEORY 

by Gerald W. Englert 
Lewis Research Center 

SUMMARY 

The early stages of a large number of avalanche formations in helium were deter- 
mined by means of computer simulation of the microscopic details. The exact trajectory 
of each electron was followed as it experienced various types of collisions with the back- 
ground atoms and as it was acted upon by the electric field in the time between collisions. 
Results of such a lengthy procedure show that the electron transport coefficients of av- 
alanches in helium relax to local equilibrium values in approximately one -hundredth of 
the formation time of a fully developed avalanche. This enables the use of a random walk 
concept for avalanche study in which the step sizes and probabilities are based on elec- 
tron transport coefficients; these coefficients depending only on the background density 
and the local electric field strengths within the avalanche. This repetitious procedure 
permits a self-consistent avalanche solution which accounts for collective long range 
Coulomb interactions as well as short range elastic and inelastic collisions with a rea- 
sonable expenditure of computer time. 

The electric field patterns initially formed inside the avalanches were close to those 
obtained from theory based on constant transport coefficients. A region of maximum and 
one of minimum induced electrostatic potential cp were located on the axis, and within 
the volume covered by the electron swarm. As formation time continues, the region of 
minimum (p moves to slightly higher radii. The electric field in the axial distance be- 
tween the two extrema builds up and diminishes in a somewhat erratic manner. 

In the intermediate time periods the avalanche growth is slightly retarded by the high 

concentration of ions in the tail which oppose the external electric field. Eventually, at a 
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value of the time-pressure parameter tp slightly greater than 3.9x10 sec-N/m 
(2.9 psec torr) the formation of ions and electrons in the regions of high field strength 
more than offset this effect, causing a very abrupt’ increase in avalanche growth. It is 
assumed that the internal electric field is axisymmetric. The random walk calculations 
were not carried much beyond the onset of the rapid growth period since assuming axi- 
symmetry is not justified there. 



INTRODUCTION 


The subject of this report is the formation of an electron avalanche, in an electric 
field, independent of wall effects (secondary ionization). This process is initially quite 
orderly, starting with one electron and continuing with an exponential growth of free elec- 
trons as a result of ionizing collisions between these electrons and the background gas. 
Such avalanches are descriptive of those which, for example, form at discrete locations 
throughout the volume of a gas, the gas undergoing initial ionization to form a plasma. 

If the formation of any single avalanche is continued under suitable conditions until suffi- 
ciently strong internal electric fields are built up, a local needle -like ionization pattern 
(streamer) is abruptly formed (ref. 1). This can ultimately result in a turbulent spark- 
ing process. This latter process is beyond the range of application of the techniques em- 
ployed in the present investigation. A fully developed avalanche will be defined herein as 
one which has grown to the streamer formation stage. 

Changes occur in avalanche formations in such short time periods that even the ex- 
ternal characteristics are difficult to measure. Theoretical solution is difficult since the 
descriptive differential equations are very nonlinear due to the extreme spatial and time- 
wise changes of the internal electric fields. A means of investigating various theoretical 
models of internal structure which allows incorporation of considerable physical detail 
could be quite useful. 

Random walk theory was employed in the computer simulation of the microscopic de- 
tails of electron motion in references 2 and 3. This theory includes both elastic and in- 
elastic short range collisions between the electrons and background neutrals , and permits 
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an exact calculation of the electron trajectory between such events. Since more than 10 
ionizations usually occur in the development of an avalanche up to the streamer formation 
stage, computer time expenditure becomes a serious problem when each electron is 
walked in detail. 

Where a large number of microscopic events occurs within small macroscopic volume 
elements and time increments, the step sizes and probabilities used in the random walk 
(RW) can be based on local averages. This in turn greatly facilitates the calculation pro- 
cedure. It was demonstrated in references 4 and 5 that such a technique can be used to 
simulate a pertinent class of parabolic nonlinear differential equations. This type of RW 
was suggested in reference 6 as a means of solution of the electron transport equation. 
The step sizes and probabilities are determined in terms of the coefficients of the trans- 
port equation which, in turn, are dependent only on the local electric field and back- 
ground density. 

In avalanche formations, the density of charged particles can reach values such that 
long range Coulomb interactions (thus internal, i.e. , induced, electric fields) must be 
considered as well as the short range interactions between the free electrons and 
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background gas particles. Application of the RW technique of references 4 and 5 to the 
electron avalanche in a manner which preserves microscopic influence, yields a self- 
consistent internal electric field, and yet enables a solution with a reasonable amount of 
computer time was undertaken. Numerical results were obtained for avalanche forma- 
tions in helium at a ratio of external electric field to background pressure of 22. 5 V-m/N 
(30 V/cm-torr). 


ANALYSIS 

Microscopic Random Walk 

In principle, the electron motion in an avalanche can be precisely simulated by ran- 
dom walks. The trajectory of each electron is followed as it undergoes various types of 
collisions and as it is acted upon by the electric field in the time between collisions. The 
details of each collision are based on random selection of values from sets of numbers 
distributed in accordance with cross section data. This is essentially the procedure of 
references 2 and 3. In these references, a constant number (representative sample) of 
electrons were followed to study a process having a uniform density of electrons. In the 
present study, however, an additional electron (therefore, walk) originates at each ion- 
izing collision, and spatial effects are of prime concern. 

Such a RW procedure provides an excellent means for computer simulation of micro- 
scopic details of avalanche formation. It requires a large amount of computer time, how- 
ever, due to the exponential buildup of electrons. This method is feasible in the early 
stages of avalanche formation (N < 10^ electrons). With this number of electrons and/or 
ions the internal electric field is negligible and only interactions of the electrons with the 
background atoms need be considered. 

This procedure was applied to the formation of a large number of electron av- 
alanches in helium. The location of the electrons and ions, at various times from the 
start of the buildup, is shown in the computer plots of figure 1. All charged particle lo- 
cations are shown here as if looking through a side view of the avalanche. That is, the 
charged particle locations are projected along lines perpendicular to the paper, which 
represents a plane through the x and z axes of the coordinate system. The origin of 
the coordinate system is the starting point (t = 0) of the trajectory of the first electron. 
The external electric field is in the negative z direction. 

The electrons essentially move in a spherical swarm. The ion pattern is that ex - 
pected from cloud chamber photographs (ref. 7). 

These avalanches were started with a number of electrons N Q considerably greater 
than one to get a higher density of points and a better definition of shape during the low 
formation times shown. The plots can also be considered to represent a superposition 
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Transverse distance, xp, N/m 


(a) Avalanche for- (b) Avalanche formation times (c) Avalanche formation time tp ■ 9. 1x10 5 
mation time tp ■ tp = 4. 5x10” 5 sec-N/m 2 sec -N/m ; number of superimposed ava- 

1,3x10”-’ sec-N/m 2 ; number of superimposed lanches N 0 = 11; number of ionizations 

number of super- avalanches N Q = 145; num- N = 1739. 

imposed avalanches ber of ionizations N ■ 1311. 

N 0 = 300; number of 
ionizations N ■ 301. 


Figure 1. - Particle locations obtained by use of microscopic RW procedure at three avalanche forma- 
tion times. Electric field parameter E/p • 22. 5 V-m/N; initial energy of first electron of N 0 ava- 
lanches is uniformly distributed from 0 to 12 eV. 


of N q avalanches, each starting with one electron, thus giving a representative average 

of formation patterns when the number of charged particles per avalanche is low. 

The electron transport coefficients were also determined in these computations, and 

are given in figures 2 and 3. Relaxation to terminal steady -state conditions were found 

to occur in short generalized time periods, tp, close to those found in reference 2, which 

is about one -hundredth of the formation time of a fully developed avalanche. Data scatter 

is due to the sample size. The number N of superimposed avalanches was selected so 
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that there was a total of about 10 electrons when formation time reached relaxation time. 

Approximately 10 minutes of computer time is required to determine formations hav- 
ing the order of 10 electrons. Avalanche structure on the order of 10 sec-N/m range 
is of main interest. Then the number of electrons may reach the order of 10*° electrons. 
Use of the detailed RW would thus require far excessive amounts of computer time during 
such long time periods. This, however, is where the use of equilibrium values of trans- 
port coefficients is justified. The electron motion can then be described by the electron 
transport equation (appendix B). This equation describes the electron motion by diffusion 
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Initial energy of 
starting electron. 
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Figure 2. - Relaxation of electron drift, diffusion coefficient, and average energy. field parameter. 

Electric field parameter E/p ■ 22. 5 V-m/N. 



and drift expressions while accounting for ionization events by a source term. 

The amount of avalanche buildup by the time that relaxation time is reached consti- 
tutes a very small fraction of the electron population in the later main time range of in- 
terest. The electron transport equation can thus be assumed to apply to avalanche forma- 
tions from the start with negligible consequence. 


Macroscopic Random Walk 

Random walk theory . - The nonlinear parabolic form of the electron transport equa- 
tion is similar to that of equations solved by the RW-on-a-grid concept in references 4 
and 5. In this RW procedure the required step sizes and probabilities are expressed in 
terms of the coefficients of the partial differential equation to be solved. 

The general equation describing a RW on a grid is (ref. 4) 



The symbols are defined in appendix A. This equation, based on the law of compound 
probability, is essentially a law of conservation of particles as they move to and from 
neighboring points on a three-dimensional lattice (fig. 4). The components of the step 
sizes AXj are the grid spaces. The terms Pj 1 " and Pr are the probabilities of taking 
steps in the positive and negative directions along the coordinate. It is assumed that 
each of the three components of the steps contribute to the particle movement during each 
time increment At. Then 


p+ + p: = i j = 1,2,3 (2) 

which corresponds to steps across diagonals of the lattice. 

Let Xj where j = 1,2,3 be Cartesian coordinates. Then W, which is the prob- 
ability density of particles being at a certain grid location at time t, can simply be set 
equal to the density of electrons n g . 

The components of step size AXj -are identifiable with second moments (root mean 
squared averages) of the components of the microscopic distance Ax^ between collisions. 
Letting a bar denote average, 




Figure 4. - Random walk on three-dimensional rectangular grid. Division 
of space into larger cylindrical zones for tallying purposes (appendix 0 
is shown by dashed lines. Radial distance r is measured from axis of 
avalanche. Axial distance z is measured from initial Ixation of first 
electron. 


Then using the definition of diffusion coefficient, 



2 r 


it follows that 


AX j = 3 = 1,2,3 (4) 

The mean free time between collisions r is equivalent to At. For helium, r can be ap- 
proximated by (ref. 2 ) 


1 _ 

(nQ a v) 


(5) 
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This is because the product of electron velocity v with Q a is very nearly constant over 
the range of interest herein. Here Q a is the total absorption cross section for electrons 
colliding with neutral background (helium) atoms of density n. Setting v equal to 
y 2e/m, where e is the mean electron energy of reference 2 for the macroscopic RW, 
gives the very small variation of r with E/p shown in figure 5. 



Figure 5. - Characteristic quantities of use in 
macroscopic RW. (Data points were calculated 
from transport coefficients of ref. 2 ; lines 
denote empirical curve fits used in the com- 
puter program. I 


The bias P? 


1 


P." can be used to describe the influence of a preferred direction of 


(drift) velocity, which involves a first moment of Ax^ 




V ^) 2 


j = 1,2,3 


( 6 ) 


Diffusion in the presence of an electric field has a nearly isotropic coefficient over 
the range of E/p of most concern in this investigation (see table I of ref. 2). It will be 
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assumed, therefore, that mean free time is constant and that step size is isotropic and 
only weakly dependent on space, although it can be strongly time dependent. Using equa- 
tions (4) and (6) in (1) and identifying W with n g then yields 


3n 


dt 


- = - V • n v, + D V 2 n 
e d e 


(7) 


Consider ionizing events to occur during the RW with a probability per collision 
P ion . The rate of ionizing events per unit volume is then P i on n e /T. This is equivalent 
to a T v d n e where avp is the Townsend first ionization coefficient. Thus 

P ion = “r v d T 


Equation (7), extended to include this source term, becomes 

V • n e v d + D V 2 n e + Qir v d n e (8) 

This is equivalent to the electron transport equation (see appendix B) within the spatial 
restriction on step size (and thus on D for constant r) which will be appraised further 
in the section entitled Results of Macroscopic RW. The transport equation- can be written 
in terms independent of background density. To follow convention, however, pressure p 
referenced to 0° C will be used as 




9n v , . 
£ d,] 

0(X jP ) 



+ 7 Ve 


(9) 


Step size and probabilities are also plotted in figure 5. These values are based on 
the transport coefficients of reference 2. 

Execution of the random walks . - In the execution of each step of the RW-on-a-grid 
procedure four numbers are drawn at random from a set of numbers uniformly distrib- 
uted over the interval 0 to 1. Such sets of random numbers are available in most com- 
puter libraries. Three of these numbers are used to determine step direction. For ex- 
ample, if the jth random number is less than the corresponding Pj 1- , the component of 
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the step along the jth coordinate is considered positive. Here, for convenience, Pt 
can be written as 



by combining equations (2) and (6). The fourth random number is compared with P ion ; 
and, if it is less than P ion , an ionization event is assumed to occur. 

Only the electrons have appreciable diffusive motion or cause ionization, and are, 
therefore, random walked. It is sufficient to include merely the drift motion of the ions 
in the calculations. 

At various times, t^, along a macroscopic time scale, the locations of the electrons 
and ions are recorded for plotting (as in fig. 1) and also tallied (similar to that in ref. 6) 
to determine their distributions in space. From the tallies, the internal electric field 
can be determined (see appendix C) and then the probabilities and step sizes from infor- 
mation as on figure 5. These quantities are determined over the physical space of im- 
mediate interest and stored in matrix form for convenient access during the walks per- 
formed in the next period ^+1 " V The appropriate (j th ) component of the total (exter - 
nal plus induced internal) electric field vector is used in the determination of the vector 
components Pj and AX^ , whereas the magnitude of the total electric field vector is used 

in the selection of the scalar P. 

ion g 

In the early formation period of the avalanche (N g < 10°) the internal electric field is 
negligible and the electron transport coefficients, thus probabilities and step sizes, are 
constants. As the formation proceeds, the permissable lengths of the time periods At^ 
become less in order to obtain a self-consistent description of collective interaction of 
charged particles with the field. As the magnitude of the internal electric field reaches 
and goes beyond that of the external field, the At^ must be kept very small (of the order 
of 10 t). 

The RW can be further facilitated by scaling step size, probabilities, and r such 
that equation (9) is generated with a smaller number of (larger) steps as shown in appen- 
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dix D. Computing time can now be reduced to about 10 of that needed for the detailed 
microscopic RW; nevertheless it still becomes excessive for avalanche formation con- 

7 

taining more than 10 electrons. With such a large number of electrons in the swarm, 

however, it is not necessary to walk each one. At the higher t samplings of about 

3 ^ 

3x10 "test” electrons were found sufficient to represent avalanche growth. Each of 

such test electrons used to continue the calculations from t^ generates a subavalanche 
over the period *k+l " V The resulting distributions are weighted for the actual num- 
ber of ions and electrons being represented. In this manner the RW process can be 
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continued through the final formation periods with total computer time expenditure of the 
o 

order of 10 minutes. 

Several possibilities exist for selecting the location of the test electron at the start 
of each subavalanche. For example, the walks of a fixed number of electrons could be 
walked from each tally zone at t = t k and the resulting tallies at l k+l weighted accord- 
ing to the spacial distribution at t^. As an alternative, a certain fixed fraction of the 
electrons at t. could be continued from their respective locations (tally zones). In the 
present investigation the starting locations were randomly selected from their spatial 
distributions at t^. In these later two cases the weighting factor is a constant for any 
*k+l - t k time interval. The procedure for random selection of values from an arbitrary 
distribution when a uniformly distributed set of random numbers is provided by the com- 
puter is discussed, for example, in reference 8. 

To facilitate calculations and to reduce computer storage requirements, it was as- 
sumed that the internal electric field of the avalanche is axisymmetric, thus cylindrical 
coordinates (r, 0,z) were used for storage (appendix C). In the tallies for the distribu- 
tions of electrons and ions each element of a matrix represents the number of electrons 
or ions within an annular zone covering distances ±Ar/2 and ±Az/2 about a point (r, z). 

A flow diagram which gives a word description of the main features of the computing 
procedure is given in figure 6. 


RESULTS OF THE MACROSCOPIC RANDOM WALK 

Only one avalanche will be presented for the RW-on-a-grid type calculation where 

long formation times, extended to the streamer stage is of interest. This avalanche, 

based on transport coefficients, has an initial growth equal to that obtained by an aver- 
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aging of 10 formations calculated by means of the microscopic RW method. Diming the 
remaining formation time, where many electrons are in the statistical sample, the mac- 
roscopic RW continues to generate the structure close to that of an average avalanche. 

A plot of the ion locations in such a fully developed avalanche, at an external E/p 
of 22. 5 V-m/N is presented in the first part of figure 7. Here the avalanche profile is 
obtained by plotting z against radial distance r, instead of x as in figure 1. This 
gives a better definition of the avalanche profile as it makes more points visible in the 
outer boundary region due to the scale factor r in the following equation: 




0 20 40 


rp = 4 V( Dp)(zp)/v d 

J Axial distance ±4yiDp)(tp) 

from v^tp (denoted by arrow) 



02040 02040 02040 


Radial distance, rp, N/m 


(a) Fully developed 
ion formation at 
time parameter 
tp- 3.9MCT 4 
sec-N/m . 


(b) Electron swarm 
at time parameter 
tp = 3. 4xl0' 4 sec- 
N/rrr ; number of 
subavalanches 
plotted N sub = 
2788. 


(c) Electron swarm 
at time parameter 
tp = 3. 68xl0“ 4 sec- 
l\l/m^; number of 
subavalanches 
plotted N sub = 
7827. 


(d) Electron swarm 
at time parameter 
tp- 3.87X10' 4 sec 
l\l/m ; number of 
subavalanches 
plotted N sub = 
4513. 


Figure 7. - Particle locations obtained by macroscopic RW procedure. Shown 
first is a fully developed ion formation; next are shown electron swarms at 
three values of time parameter. External electric field parameter E/p ■ 

22. 5 V-m/N. 


Twice the diffusion radius, 2r^ = 4 |/Dz/v^, shown by the dashed line fairly well defines 
the avalanche boundary. It defines the boundary within which 1 - 1/e^ (or 98 percent) of 
the electrons are located. This can be easily obtained by use of equation (B4) of appen- 
dix B. 

The corresponding electron swarm is shown at three time periods. It is nearly 
spherical at the first time period since the internal electric field is then quite small. It 
has a radius close to 2r ^ and its center is close to the axial location z = v^t. By the 
end of the second time interval the internal field reaches values close to half that of the 
external field and now definitely influences the electron pattern. Local minimum and 
maximum values of internal E occur at zp values of approximately 77 and 100 N-m, 
respectively. The concentration of electrons is strongly correlated to such field points 
since avp, and thus the formation of new electrons and ions, is exponentially dependent 
on field strength. 

4 2 

As the avalanche buildup continues to tp slightly greater than 3. 6x10 sec-N/m , 
the electric field strength near the axis and in the vicinity of zp = v^tp increases to mag- 
nitudes beyond the external field strength. The last plot on figure 7 shows that there is a 
corresponding large concentration of electrons (and, therefore, ions) at this location. In 
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fact, the density is so high at tp = 3. 9xl0' 4 sec-N/m 2 that only a small sample of elec- 
trons can be plotted without a glare in this zone occurring on the plotting screen and ob- 
scuring details. 

Contour plots of the ion and electron densities are shown in figure 8 at a formation 
-4 2 

time of 3.25x10 sec-N/m . The RW results (fig. 8(a)) compare quite well with con- 
stant coefficient theory (fig. 8(b)) at this tp where the induced electric field is much 
less than the external field. The diffusion radius profile (2r D ) is also shown on the RW 
ion plot. The exponential falloff of ion, as well as electric, density from the peak re- 
gions near z = v^t and r = 0 is apparent. 

Marginal distributions, in which integration is over all of the dimensions except the 
one being studied, are used to present the effects of the internal electric field on density 
distribution. These often allow a space saving in plotting while preserving the essential 
features of main interest. Such plots are shown in figures 9 and 10 for both charge spe- 
cies at various formation times. The marginal density distributions (n. and n. ) for 

j,r j,z' 



(a) Time parameter tp • 2. 13x10 4 sec-N/m 2 ; (b) Time parameter tp ■ 3. 25x10 4 sec-N/m 2 ; 

effective Townsend ionization coefficient effective Townsend ionization coefficient 

(Off/pieff ■ ay/p = 0. 238 ml N. ^P> e ff ’ °\ f P = 236 m,N - 



0 4 8 12 15 18 21 2d 0 4 8 12 15 18 21 24 

Radial distance, rp, N/m 


(c) Time parameter tp = 3. 68x10 " 4 sec-N/m 2 ; (d) Time parameter tp = 3. 89xl0' 4 sec-N/m 2 ; 

effective Townsend ionization coefficient effective Townsend ionization coefficient 

(Of/p) eff = V P = 0, 234 m/N. (or/pleff = Ojlp - 0. 229 m/N. 

Figure 9. - Marginal density distribution in radial direction at various values of time parameter. 
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(c) Time parameter tp ■ 3.68x10'“ sec-N/m 2 ; 
effective Townsend ionization coefficient 
(Of/plgjf = dp/p = 0. 234 m/N. 


M 



(d) Time parameter tp = 3. 89xl0' 4 sec-Nfm 2 ; 
effective Townsend ionization coefficient 
(Of/p) eff = Of/p = 0. 229 m/N. 



(f) Time parameter tp = 3. 908xl0" 4 sec-N/m 2 . 


Figure 10. - Concluded. 


the j th charge species in the r and z directions, respectively, are normalized by 
dividing by exp( acpV^t) yielding 

r °° r2v 

F j,r =n j r exp(-cn r v d t) = exp^a^v^) J n.(r,e,z)r de dz (12) 

and 

/~<x> r2n 

F i z = n j z expC-Q^Vjjt) s expi-a^v^t) n.(r , e, z)r de dr (13) 

Keeping the scale factor r with the integration over 9 in equation (12) makes the plots of 

F- against r more manageable as it reduces the otherwise exceptionally high values 

near r = 0. Theoretical values of n. „ and n. _ are given in equations (B7) to (B10). 

J ) * J) “ 

When the internal electric field is negligible, so that the transport coefficients are 
constant, theory agrees well with the RW results, as expected. The ion distribution 
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peaks at a very slightly lower radius than the electrons (fig. 9). A difference between 
the distribution of the two charge species is more noticeable in the axial direction due to 
the electron drift velocity. 

-4 / 2 

At about a formation time of 3. 2x10 sec-N/m the difference between the RW and 
the approximate theory becomes noticeable; especially, for the z component of the dis- 
tributions (fig. 10(b)). The theory is shown only for the r component for the next two 
time periods selected for plotting; beyond which it is inapplicable. The lines on the sub- 
sequent plots are merely fairings through the RW data. 

-4 2 

Near tp of 3.6x10 sec-N/m the z distribution curve starts to neck down at a z 
location somewhat less than v^t, and bulges out on either side of this z. This is evident 
in figure 10(c). Note that the computer plots of figure 7 show a tendency for the electron 
swarm to separate into two zones. The corresponding radial distribution (see fig. 9(d)) 
now shows an increased concentration of both charge species near the axes as well as a 
spreading of the outer parts of the curve to radial distances beyond the theoretical curve. 
As time continues, these trends become more pronounced with a sharp peak forming near 
the axis and near v^t in the radial and axial distributions, respectively. This zone is 
apparent in the last computer plot of figure 7. By this time the induced electric field has 
exceeded the magnitude of the external field in this region. 

Contour plots of the induced electrostatic potential cp are shown in figure 11. The 
lines of constant cp form sets of approximately circular contours about a region of max- 
imum <p located at a z less than v d t and a region of minimum cp at a z greater 
than v^t. The cp - 0 line is located near a z equal to v^t during the early formation 
period. 

The internal electric field in the avalanche head is generally predicted to be directed 
inward toward the center (pp. 78 and 79 of ref. 1) forcing the electrons outward; whereas 
the internal field in the trailing zone, being rich in ions, is predicted to be outward tend- 
ing to force the electrons inward. This general pattern is observed in the present study 
only during the early formation period where the extrema of the field strength are both 
located on the avalanche axis. 

Comparison of figures 11(a) and (b) show again that the constant coefficient theory 
(eqs. (B12) and (B13)) agrees quite well with RW results at low tp where internal fields 
are much less than the external value. The contours obtained by both methods look quite 
alike over much of the field; however, it is the local extrema regions which are decisive 
in the latter avalanche behavior. A plot of the values of cp extrema against time is 
shown in figure 12. 

4 2 

As formation time continues beyond 3. 6x10 sec-N/m , the region of minimum cp 
moves to slightly higher radii at times (fig. 11(c) and (d)). Comparison with figure 6 
shows that the two extrema regions and the distance between them is essentially the space 
covered by the electron swarm. 
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(e) RW results at formation time 
tp - 3.907xl(T 4 sec-N/m^. 

Figure 11. ■ Concluded. 



Figure 12. - Buildup of maximum and minimum induced 
potentials with time. 


The contour lines about and in the near vicinity of the regions of extrema often have 
two or more inflection points oriented such that they indicate enhanced gradients of <p in 
the radial direction. During the latter stages of avalanche formation, the radial electric 
field component is found to grow more rapidly than the axial component. A plot of E r /p 
against axial distance in the regions of extrema is shown in figure 13 at two radial loca- 
tions near the axis and at various formation times. 

At low tp the radial electric field has approximately a sinusoidal variation as zp 
increases over the span of the electron swarm. The variation of E r with z at the inner 
radial location of figure 13 is more erratic, and indicates a tendency of local charge 
buildup and relaxation at z values between the two initial extrema (see also fig. 11(e)). 
This may be an indication of instability caused by the exponential dependence of a^/p on 
E/p. This causes a greatly increased rate of electron and ion formation at high E/p lo- 
cations. The restoring, or stabilizing influence is primarily the drift of electrons in re- 
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Negative radial component of electric field strength, -E^p, V-m/N 
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■ 3. 72xl0" 4 sec-N/m 2 . 


Radial distance, 

rp, 

Nlm 

1.10 

3.33 



Axial distance ±4 V< Dpi(tp) 
from vJp (approximate 
axial distance covered by 
electron swarm) 





Axial distance, zp, N/m 
(g) Formation time tp = 3. 90xl0" 4 sec-N/m 2 
Figure 13. - Buildup of radial components of internal electric field extrema. 
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sponse to the local electric field. At a formation time of 3. 91xl0" 4 sec-N/m 2 the field 

at the inner radius increases so abruptly that the assumption of axial symmetry can no 

longer be justified and the calculations were terminated. This abrupt change occurs at a 

03 

formation of about e ions which somewhat surpasses the estimate of e u for break- 
down in reference 1 (p. 131). 

Another indicator of the influence of the induced internal electric field is the effective 
Townsend ionization coefficient, defined here as 



A retardation of avalanche growth below that calculated with a constant a^/p equal to 
( Q! q'/p)( ;= o predicted in references 9 and 10. Such an effect is observed in the present 
results, as shown in figure 14(a), over about the latter one-third of the formation time. 



fb! Size dependence. 

Figure 14. - Influence of induced internal electric field on effective Townsend ionization 
coefficient. 
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This is due primarily to the ions lagging behind the electrons and causing an internal 

_4 

electric field which opposes the external field. At a tp slightly beyond 3. 87x10 sec- 
o 

N/m however, the high field strengths discussed previously cause local increases in 
<* T /p which more than offset the retardation process and (aip/p)^ sharply increases. 

A plot of ( Q ^/p) e ff against N is given in figure 14(b). The general trends are quite 
consistent with the experimentally obtained results of reference 11 for Ng. 

Use of the (a^/p)^ of figure 14 in place of (c^/p^g in equations (B8) and (BIO) 
gave little change in the theoretical ion marginal distributions of figures 9 and 10. Use 
of ( a 'p/p) e ff in the electrostatic potential calculation of equations (Bll) and (B12), how- 
ever, noticeably improved agreement between approximate theory and the RW results in 
the intermediate time periods. 

The assumption that step size has only a weak spatial dependence, used in the der- 
ivation of equation (7) can now be further appraised. The ratio of step size AX. to the 
step size at z = 0 is plotted in figure 15 against zp at a low value of rp and a high 
value of tp where the internal field, and thus transport coefficients, undergo extreme 
variations. This is compared to the dependence of other quantities entering the macro- 
scopic RW calculations: the z component of bias, P* - P“, and the probability of ion- 
ization P ion . Step size is shown to be quite insensitive to distance, especially compared 
to the large changes experienced in the other two RW quantities. 



Axial distance, zp, N/m 

Figure 15. - Spatial variation of step size AXj, z component 
of bias P* - PjJ, and probability of ionization Pj on . Radial 
distance from avalanche axis rp = 1.1 N/m; formation time 
tp - 3. 8 sec-N/m Z ; external electric field strength E/p = 

22. 5 V-m/N. 
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CONCLUSIONS 


The present random walk (RW) study of avalanche formation in helium leads to the 
following conclusions: 

1. The results of random walks utilizing simulation of the electron trajectories in an 
avalanche in minute detail show that the electron transport coefficients very quickly relax 
to local equilibrium values. This enables the use of an alternate RW-on-a-grid concept 
in which the step sizes and probabilities are based on equilibrium transport coefficients. 
Thus they are functions only of the local electric field and background density of neutrals. 
The concomitant reduction in computer time expenditure makes calculations of large av- 
alanche formations feasible. By only a simple repetitious procedure a self-consistent in- 
ternal electric field solution is obtained. Collective long range Coulomb interactions as 
well as short range elastic and inelastic collisions are accounted for. 

2. Results obtained by this method at an E/p of 22. 5 V-m/N show that the internal 
electrostatic potential cp pattern in its early stage of formation has minumum and max- 
imum values of potential on the avalanche axis. The axial distance over which large cp 
gradients appear is essentially that spanned by the electron swarm. As avalanche growth 
continues, the minimum <p region moves to slightly higher radii. Within the separation 
distance between the two extrema, several other large ^-gradient zones form and die 
out in a manner which appears to be random in time. At a generalized formation time tp 
slightly greater than 3. 87x10 sec-N/m, the internal electric field abruptly increases 
far beyond the magnitude of the external field. 

3. The tail region behind the electron swarm, has a positive charge which opposes 
the external field. This tends to reduce the effective ionization coefficient, which is an 
indicator of the avalanche growth. At the time of the abrupt increase of the magnitudes 
of the internal field strength, the electron and ion formation rate in the region of the (p 
extrema increases to the extent that it offsets the growth retardation influence in the tail 
region. The effective ionization coefficient then sharply increases. 

4. Theory, based on constant transport coefficients, gives good results only in the 
early formation periods of the avalanche where internal electric fields are much less than 
the external fields. Here, electron and ion density distributions, as well as the general 
contours of the internal electrostatic potential are quite satisfactory. Some noticeable 
improvement in the theoretical <p contours is obtained by replacing Townsend’s first 
ionization coefficient based on the external field by the effective ionization coefficient for 
the formation time of interest. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 4, 1973, 

503-10. 
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APPENDIX A 


SYMBOLS 

A scale factor 

a combination of variables defined by eq. (C7) 

b combination of variables defined by eq. (C8) 

D coefficient of diffusion 

E electric field strength 

£ elliptic integral of second kind 

e electrostatic unit of charge 

F. r normalized marginal density distribution in radial direction 

F. normalized marginal density distribution in axial direction 

J> z 

& elliptic integral of first kind 

X complete elliptic integral of first kind 

N number 

n particle density (with no subscripts this refers to neutral background particles) 
n. marginal density distribution of species in radial direction 

J 9 1 

n. marginal density distribution of jth species in axial direction 
J > z 

P ion probability that encounter between electron and neutral atom causes ionization 
pt probability that j** 1 component of step will be in positive direction 

Pr probability that jth component of step will be in negative direction 

p pressure 

total absorption cross section of electrons colliding with background neutral atoms 
R radial distance in spherical coordinates 

r radial distance, in cylindrical coordinates, from axis of avalanche 

Tp diffusion radius 

S source density 

t time 

v velocity 
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V 

W 

x i 

X j 


a, 


T 

y 

6 

€ o 

7 r 

V 

0 

e 

£ 

p 

T 

* 


mean velocity obtained by integrating over distribution of total velocity 
probability density 

distance in jth direction used in macroscopic descriptions 
distance in jth direction used in microscopic descriptions 


x, y,z Cartesian avalanche coordinates; x and y are in a plane perpendicular to z; 
z is axial distance along avalanche axis 


Townsend first ionization coefficient 
combination of variables defined by eq. (C9) 
combination of variables defined by eq. (C 13) 
small finite increment 
capacitivity of vacuum 

mean electron energy obtained by integrating over distribution of random energies 
dimensionless distance in avalanche space 

angle between field and source vectors in cylindrical coordinates 
angle 

radius ratio 
charge density 
mean free time 
induced electrostatic potential 
supplement of angle 0 
Subscripts : 
a absorption 
diffusion 
drift 
electron 
ion 

ionization event 

! 

I 

initial 
radial 
axial 


D 

d 

e 

i 

ion 

o 

r 

z 
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APPENDIX B 


DERIVATION OF MACROSCOPIC DESCRIPTIVE EQUATIONS 

For the range of electric field strength of interest here, a diffusion coefficient for an 
isotropic medium can be used (see table I of ref. 2). The flux of electrons about a mean 
(drift) velocity v^ due to diffusion is then given by Fick's law as (ref. 12) 

n e^- f d ) = - DVn e 

Substituting this into the continuity equation, 

3n ■» — 

— — + V • n v = S(r) 

at e 

where S(r) = accounts for sources, gives the electron transport equation 

Sn_ 

— = “ v * n e v d + V ‘ (D Vn e ) + Q! T v d n e (B1) 


Here is the Townsend first ionization coefficient and o^v^^ is the increase of the 
number of electrons per unit time per unit volume due to ionization events. Only singly 
charged ions are considered. 

The ions are cold and have a negligible random velocity. Their transport equation 
is simply 

Sn. 

"i T d,i + “r T d n e < B2 


Equations (Bl) and (B2) are coupled through Poisson's equation, 


2 -e(n, - n ) 

v <p = 


(B3) 


and the dependence of the transport coefficients on the electric field vector which is equal 
to - V(p plus the external electric field. 
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Density Distributions 

For constant coefficients equations (Bl) and (B2) can be solved immediately. The 
electron density is 


n e 


(47rDt) 


3/2 


exp 


■w - 


r 2 + 


( z - y) : 


4Dt 


(B4) 


in cylindrical coordinates with azimuthal symmetry. Here the avalanche is assumed to 
start from one electron. The ion density can then be obtained from 


n. - 


a T v d 


l‘ n ‘ 


dt’ 


(B5) 


if it assumed that the ion drift is negligible. Using equation (B4) in equation (B5) gives 


n. = 


¥d 


t- 

J 0 (477Dt-) 3/2 


Vd v - 


r 2 + 


( z - V') 


A2' 


4Dt' 


dt’ 


(B6) 


Equations (B4) and (B6) are used in the avalanche studies of references 9 and 13. 

Marginal distributions are often used to reduce the number of variables required for 
plotting results. Integrating densities, n e and n., over r and 6 gives the marginal 
distributions, n _ and n. . in z as follows. By definition 

C j Z 1 j Z 7 


r°° r2ir 

n e,z=7 0 J 0 ”e rd « dr 


Substituting equation (B4) for n g yields 


n. 


2n 


6,Z (4M) 3 /* 


exp 


“tV 


(z - T d , 


4Dt 


l]f 

- 


r exp )dr 

\ 4Dt/ 


yiffDt 


exp 


rrV 


( z - v)’ 


4Dt 


(B7) 
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Using this result in equation (B5) yields the following expression for ions: 


Ctr 


n. 


l,z 


f i 

V* J 0 yiot 7 


exp 


A2' 


W 


( z - V) 


4Dt* 


dt’ 


(B8) 


In like manner 


where 


r°° 

"e,r-J 0 l V d « dz 


27 T yjptr 
(47rDt) 3/2 



(z - v d t) 

V4Dt 


Finally 


n e r = — exp (a r v d t - — ) 1 + erf 
’ 4Dt V Td 4Dt/[ 

For formation times of interest erf(v d t/y4Dt) « 1, and 


n 


e,r 


r 

= — exp 
2Dt 



Using this result in equation (B5) gives 



Equations (B7) to (BIO) can be normalized by dividing by exp(aj,v d t). 


(B9) 


(BIO) 
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Potential Distributions 


Use of the assumption of constant transport coefficients in the preceding derivations 
results in expressions which show that the electrons then diffuse as a spherically sym- 
metric swarm, the center of which is at z = v^t on the avalanche axis of symmetry. 
This simplification in turn permits determination of the electrostatic potential cp Q due 
to the electron distribution in space. 

The electrostatic potential for electrons can be obtained from 


<P e ( r > t) 




d 3 r' 


where integration is over the whole field of the electrons. In like manner the potential 
for ions is 



Using equation (B5) and singly charged particles 



so that 


<p.(r,t) = -a r v d 



<P e (r,t’)dt' 


(Bll) 


The total potential <p is then the sum of cp e and <p^. 

Let R denote radial distance in the spherical (avalanche space) system. The well- 
known potential equation for such a charge distribution is 


<P e (R,t) = 


/ R p°° 

R’ 2 p e (R*,t)dR’ + — / RV R ’> 

€ ° Jr 


t)dR’ 
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Using equation (B4) gives 


<pJR,t) - - 


e exp(oi T v d t) 
, 3 / 2 e 0 |R 



R/R 


D 


d« + 


R 



«e'« d« 


R/R. 


D 


where £ = R/R D and R Q = |^Dt. Finally 


e exp(d! T v rI t) / R \ 
cpj R,t) = -LA_ erf (A 


4?re o R 


R 


D/ 


where the method of integration by parts was used in evaluation of the first integral. 
Changing back to cylindrical coordinates in laboratory space gives 


1 


<P e (r, z,t) = 


e exp(o! T v d t) 


erf 


r 2 + 


( z - v) : 


-, 1/2 


4ne Q Vr 2 + (z - v d t) 2 

Substituting equation (B12) into equation (Bll) then yields 


4Dt 


(B12) 


<pft, z,t) = 


Q4pV d e I expCo^v^’) 


4tt€ " 


Vr 2 + (z -v d t') 5 


. erf 


.2 . L .. .,\2i l/2 


r + 


(■ - vf 


4Dt' 


dt' (B13) 


^he radial electric field E R is now given by 


E R = 


d<p _ e exp(a T v d t) d 
3R 


4tt€ R 2 d — 

o D R 

R D 


erf/A N 

\ D/ 


R 


R 


D J 


which agrees with equation (4) of reference 9. 
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APPENDIX C 


DETERMINATION OF INTERNAL ELECTRIC FIELD 
FROM TALLIES OF ELECTRONS AND IONS 

Axial symmetry is assumed so that the induced electric field has only a radial, E r , 
and an axial, E z , component. The field and source points are located at r, 0 ,z and r', 
e\ z', respectively. Using Coulombs law, and MKS units, the induced electrostatic po- 
tential for singly charged ions is 


<p(r,z) = 


e f f f n t (r’,z') -n e (r',z’ 

47r€ ° J J J Vr 2 + r' 2 - 2rr' cos 0 + 


r ’ dr ' d© dz 


(Cl) 


(z - z’) 2 


where 0 = 6 - 6' and the integration is over the whole space occupied by the avalanche. 
By differentiation it follows that 


E p (r,z) = 


e fff [ n i (r ’’ z,) " n e (r',z’)](r - r» cos 0) ( ^ 

° J J J |r 2 + r' 2 - 2rr T cos 0 + (z - z’)^] 


i3/2 


r ' dr' d0 dz (C2) 


and 


E z (r,z) - 


^ [f[_h 

4u oJJJ r r 2 r , 


[n^r', z’)- n e (r’,z')J(z - z’) 

3/2 


r’dr'd0dz (C3) 


[r 2 + r' 2 - 2rr' cos 0 + (z - z’) 2 


The integrations over 0 can be expressed in terms of elliptic integrals. As long as 

r' and z’ do not equal r and z simultaneously, equations (Cl) to (C3) can be reduced, 

for example, by using equation (5) page 154 and equation (4) page 156 of reference 14. 

Indicating the contribution of the source points to <p, E r , and E z in the neighborhood of 

the field points by 5cp, 6E„, and 6E„ and integrating over the remaining source volume 

r z 

gives 
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cp = 6<p + 


” £ o JJ fiTi 


jf(r)r' dr' dz' 


(C4) 


E r = 5E r + 


7T€, 


/ n i " n e 

br - ar» + 

/ b ^a + b 

a - b 


r' dr' dz’ 


(C5) 


and 


E z = 6 E z + 


ff ^ "I, 

f€ o J J (a - b) y'a + 


(z - z')<?(r)r’ dr' dz’ 


(C6) 


where 


+ r’ 2 + (z - z’) 2 

(C7) 

b s 2rr’ 

(C8) 


and 


r = 


V 2b 

a + b 


(C9) 


The identity 

cos 6 _ 1 a 1 

(a - b cos b f^a - b cos 0 13 (a - b cos O)^ 2 

was used in simplifying equation (C2). 

When r approaches r' and z approaches z r , slightly more cumbersome expres- 
sion result. The distance of closest approach between charged particles is assumed to 
be n e -1// ^. It is now convenient to replace 0 with its supplement, ip = n - 0. Then, 
using, for example, equations (4) and (3) of pages 154 and 156, respectively, of refer- 
ence 14, results in 
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e 


(CIO) 


6 (p = 


2?re, 


JJ [n^r^z’) - n e (r’,z')]y| jr^y, ijr’ dr’ dz' 


5E r = 


i t a JfhW'W\)l Irrf 4 r! 


2(rb - ar') 
(a - b)(a + b) 


£ 


i 4 - ? 


sin 6i \j 


^a - b cos 


>r’ dr' dz’ 


6Ez “^o / f r) 


2b(z - z’) 

(a - b)(a + b) 


£ 


- *(y, - 

b 


sin 5x}y 


>r' dr' dz' 


r / |^a - b cos 5i/J_ 


where 


+ cos 6i^) 
a + b 

Here the distance of closest approach in the \f/ direction is expressed as arc length 
r6^, or 


y = arc sin 


i 


6^ = 


rn 


1/3 


As r approaches r' and z approaches z', r approaches 1 and 


y 


-sin _1 ^/I± 


cos 6i^ 


(Cll) 


(C12) 


(C 13) 


(C 14) 
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Then 


1 


— - S(y, 1) = sin y = 

v TV \ 




cos 6^ 


and 


y, I = In tnfl. + — J 

V r [2 4 


Thus for small separation distance between field and source points 


5cp = 


5E r = 


~ jj[ HiCr’.z’) - n e (r',z')J 

~ jj [ n i (r ’’ z,) - V r ’’ z,) ] 


In tn [ — + — )r' dr’ dz’ 
\2 4 / 


In tnf ^ + — |r’ dr' dz' 
\2 4 


and 


6E z = 


JJ n- 1/3 [ ni (r', z ') - n e (r',z')]lntn^ + l' 


r’ dr’ dz’ 


(C15) 


(C 16) 


(C 17) 


In the computing procedure the ion and electron density distribution in r,z space is 
represented in matrix form. Each element of a density matrix represents the number of 
ions or electrons inside a thin or small annular source volume (zone) equal to 
2irr' Ar’ Az’ (fig. 4). The number of particles N(i,j) in each zone is obtained during the 
random walks by a tally procedure at the t^ time periods of interest. The cp, E r , and 
E are assumed constant in each of such small field zones. Using this approximation in 
equations (C4) to (C17) yields 


(p{l , m) = 




N^m) - N e (Z,m) 


r" 1 In tn Hi + -j 


+ 2 


20 

I 

j,k=l 


* [NjG,k) -N e Q,k)] 
^a + b 


*(r)l 


(C 18) 
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' ntn \2 + 4j 


20 

V’ NjQ.k) -N e (j, 

jfei 


k) 


br - ar* 
a - b 


<£(r) + r*jf(r) 


(C 19) 


and 


E z (Z,m)=— ® J N.(Z , m) - N e (Z , nOjr^n" 1 / 3 In tn (l. + l\ 

16tt 2 € I \ 2 4 / 


+ 8 


Z ’N.(j,k)-N e (j,l 

, (a - b) t/a + b 

],k=1 * 


(z-z’K(r)y (C20) 


where the prime marks on the summation signs indicate exclusion of j = l and k = m 
combination. The r, r*, z, and z' distances which enter the equation are measured to 
the midpoints of the annular zone; that is, the inner radii + Ar/2 and to the lower z 
boundary + Az/2. Limitation of each group of zones to 400 was found to provide suffi- 
ciently accurate charge distributions and yet be within the limits of a convenient amount 
of computer storage. 
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APPENDIX D 


SCALING THE MACROSCOPIC RW 


In the early time periods when the avalanche is undergoing relatively slow changes, 
the calculation procedure can be further facilitated by using a scale factor, A. Let 



-p:’= 


P ion - AP ton 


and 


t' = At 


Then using equations (4) and (6) 



and 



From the definition of P. 

ion’ 



^ion 


- 04p 


Thus the coefficients in the electron transport equation remain unchanged but the RW 
is performed with a fewer number of steps. In helium at an E/p of 22. 5 V-m/N the 

• I 

scale factor must not get much greater than 40 in order to keep PT and P| Qn less than 
unity. 
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At higher formation times the use of scale factor is of little value as it must be re- 
duced to nearly one in order to keep r* and AX.' much less than the time and spatial 
increments within which there are sizable avalanche changes. 
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applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 




